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Summary: Scaling predictions and results from self- 
consistent field calculations for bottle-brush polymers 
with a rigid backbone and flexible side chains under good 
solvent conditions are summarized and their validity and 
applicability is assessed by a comparison with Monte 
Carlo simulations of a simple lattice model. It is shown 
that under typical conditions, as they are also present 
in experiments, only a rather weak stretching of the side 
chains is realized, and then the scaling predictions based 
on the extension of the Daoud-Cotton blob picture are 
not applicable. 

Also two-component bottle brush polymers are consid- 
ered, where two types (A,B) of side chains are grafted, as- 
suming that monomers of different kind repel each other. 
In this case, variable solvent quality is allowed for, such 
that for poor solvent conditions rather dense cylinder- 
like structures result. Theories predict "Janus Cylinder" - 
type phase separation along the backbone in this case. 
The Monte Carlo simulations, using the pruned-enriched 
Rosenbluth method (PERM) then are restricted to rather 
short side chain length. Nevertheless, evidence is ob- 
tained that the phase separation between an A-rich part 
of the cylindrical molecule and a B-rich part can only 
occur locally. The correlation length of this microphase 
separation can be controlled by the solvent quality. This 
lack of a phase transition is interpreted by an analogy 
with models for ferromagnets in one space dimension. 



I. Introduction 

Flexible macromolecules can be grafted to various sub- 
strates by special endgroups.|--|Su|eh|--|'j8olymer brushes" 
find widespread applications'^ l^'l^' 1^1^ and also pose 
challenging theoretical problems, such as an understand- 
ing of the conformational statistics and resulting geomet- 
rical structure of these tethered chain molecules. Only 
this latter aspect shall be considered in the present pa- 
per, for chains grafted to a straight line or a very nar- 
row cylinder. This problem is a limiting case of "bot- 
tle brush" polymers where side chains are grafted to a 
long macromolecule that forms the backbone of the bot- 
tle brush. When this backbone chain is also a flexible 
polymer and the- grafting density is not very high, a 
"comb polymer" results, which is outside of consider- 
ation here. Also we shall not discuss the case where the 
backbone chain is very shorjtj |ao-jth*j eepformation would 
resemble a "star polymer" .'^'^''Sl>li3>liy Here we restrict 
attention to either stiff backbone chains or high graft- 
ing density of side chains at flexible backbones. In the 



lattice case stiffening of the backbone occurs due to ex- 
cluded volume interactions, and a cylindrical shape of the 
molecule as a whole results. In fact, many experiments 
have been carried out where with an appropriate chemical 
synthesis bottle brush polynirefis-witlipaq worm-like cylin- 
drical fbfiM r*pre produced li3> li3 li4> lia lial _ xhe recent 
papersli^li^'li^ contain a more detailed bibliography on 
this rapidly expanding field. 

On the theoretical side, two aspects of the confor- 
mation of bottle brush polymers where mostly dis- 
cussed: (i) conformation of a side chain when the back- 
bone caii^ pteagUjsdi ^ r*, rig'd straiffkL line or thin 
cylinderliali3'liaM|2o|.y,l2i,2;3,24, 25,y:y (ii) confor- 

mation of the wh^o— bottle-, bsusii rWjh^|^^ 
is (semi)flexible. Id. Issl, lasl, lai y . y. lai ^ y ImI M M M 
The latter problem is left out of consideration in the 
present paper. Problem (i), the stretching of the side 
chains in the radial direction in the case of sufficiently 
high grafting density, i was moatlv r-diseyssed-jn terms of 
a scaling descriptioH- ll2l> lis liSi, 123, M I25I M extending 
the Daoud-Cotton ^ "blob picture" 1^ 1^ from star 
polymerSrte rb«ttle brush polymers. If one uses the Flory 
exponent'^ 1^ v = 3/5 in the scaling relation for the 
average root mean square end-to-end distance of a side 
chain. Re cx cr(i-'^)/(i+'')iV2'^/(i+''), where a is the graft- 
ing density and N is the number of effective monomeric 
units of a side chain, one obtains Re oc a 

/4^3/4^ These 

exponents happen to be identical to those which one 
would obtain assuming that the chains attain quasi- 
two-dimensional configurations, resulting if each chain is 
confined to a disk of width a^^ . Although this lat- 
ter picture is a misconception, in experimental studies 
(e.g. [H, [1^) this hypothesis of quasi- two-dimensional 
chains is discussed as a serious possibility. Therefore we 
find it clearly necessary to first review the correct scaling 
theory based on the blob picture, and discuss in detail 
what quantities need to be recorded in order to distin- 
guish between these concepts. Thus, in the next section 
we shall give a detailed discussion of the scaling concepts 
for bottle brush polymers with rigid backbones. 

Thereafter we shall describe the Monte Carlo test 
of these predictions, that we have recently per- 
formed usittg-jtlia |pruned enriched Rosenbluth method 
(PERM). yi^'y fist After a brief description of the 
Monte Carlo methodology, we present our numerical re- 
sults and compare them to the pertinent theoretical pre- 
dictions. 

In the second part of this paper, we discuss the 
extension from one-component to two-component bot- 
tle brush polymers. Just as in a binary polymer 
blend (A,B) typically the energetically unfavorable in- 
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teraejtieni— id«S(pibed by the Flory-Huggins parameter 
^ |M|45|,15C|, [511,151^ should cause phase separation be- 
tween A-rich and B-rich domains. However, just as 
in block copolymers where A[-Gh|ains--.and B-chains are 
tethered together in a point, '^''^'^ no macroscopic 
phase separation but only "microphase separation" is 
possible: for a binary (A,B) bottle brush with a rigid 
backbone one nftajy|— expect formation of "Janus Cylin- 
der" structures. '^''^'^ This means, phase separation 
occurs such that the A-chains assemble in one half of the 
cylinder, the B chains in the other half, separated from 
the A-chains via a flat interface containing the cylinder 
axis. However, it has been argued that the long range 
order implied by such a "Janus cylinder" type struc- 
ture has a one-dimensional character, and therefore true 
long range ovdm: is destroyed by fluctuations at nonzero 
temperature.'^ Only local phase separation over a fir^t* 
correlation length along the cylinder axis may persist.'-^ 
We shall again first review the theoretical background 
on this problem, and then describe the simulation evi- 
dence. We conclude our paper by a summary and out- 
look on questions that are still open, briefly discussing 
also possible consequences on experimental work. How- 
ever, we shall not deal with the related problems of mi- 
crophase separation of a bottle brush with only one kind 
of sider-Ghains induced by deterioration of the solvest 
quality'^ 1^ or by adsorption on flat substrates. ^^iS2i>\§M 



II. Conformation of Side Chains of Bottle 
Brushes under Good Solvent Conditions: 
Theoretical Background 

The most straightforward approach to understand the 
conformations of chains in polymer brushes and star 
polymers under good solvent conditions uses the con- 
cept to partition the space available for the chains into 
compartments, called "blobs". The idea is that in each 
such region there occur only monomers of one chain, no 
monomers of any other chains occur in such a blob, and 
hence self-avoiding walk statistics holds in each blob. 
This means, if a (spherical) blob has a radius vb and 
contains n monomers, these numbers must be related via 

rB = an'' , 0.588 , n > 1 . (1) 

Here a is a length of the order of the size of an effec- 
tive monomer, and we emphasize from the start that it is 
crucial to use the correct value of the self-avoiding walk 
exponent it is provided from renormalization grcM» 
calculations'*'''' or accurate Monte Carlo simulations. 
If one would ignore the small difference between v and 
the Flory estimate 3/5, one would already miss an impor- 
tant distinction between twot-different scaling regimes for 
a brush on a flat substrate, An almost trivial condi- 
tion of this "blobology" is that each of the N effective 
monomeric units of a chain must belong to some blob. 



So we have 

N = nnB , (2) 

where ub is the number of blobs belonging to one par- 
ticular chain. 

Finally we note that the space available to the chains 
must be densely filled with blobs. It then renaains to 
discuss which factors control the blob size tb-^^ The 
simplest case is a polymer brush on a fiat substrate (un- 
der good solvent conditions, as assumed here through- 
out): if we neglect, for simplicity, any local fiuctuations 
in the grafting density, the distance between grafting 
sites simply is given by a~^/'^. Putting tb = cr^^/^ in 
Equation we find n — {aa?)~^/'^" , i.e. each chain 
is a string of ub — N/n — N{aa^y^'^'^ blobs. Accord- 
ing to the simple-Hiipded description of polymer brushes 
due to Alexander,'^ this string simply is arranged like a 
one-dimensional cigar, and one would conclude that the 
height of a fiat brush is 

h = a-'/^UB = iVa(aa2)i/(2-)-i/2 . (3) 

The free end of the chain is in the last blob and hence the 
end-to-end distAnceJl^jt^-hr-in this "Alexander picture" 
of polymer brushes. '^'^'^''^ However, a more detailed 
theory of polymer brushes, such as the, self-consistent 
field theory in the strong stretching limitj'S^I'Issl.lssI yields 
a somewhat different behavior: the end monomer is not 
localized at the outer edge of the brush, but rather can be 
located anywhere in the brush, according to a broad dis- 
tribution; also the monomer density in polymer brushes 
at flat substrates is not constant up to the brush height 
h, but rather decreases according to a parabolic profile. 
So even for a polymer brush at a flat substrate already 
a description in terms of non- uniform blob sizes, that 
increase with increasing distance z from the substrate, 
is required.'^ However, in the following we shall disre- 
gard all these caveats about the Alexander picture for 
flat brushes, and consider its generalization to the bottle 
brush geometry, where polymer chains are tethered to a 
line rather than a flat surface. Then we have to partition 
space into blobs of nonuniform size and shape in order to 
respect the cylindrical geometry (Figure [1]). 

In the discussion of brushes in cylindrical geometry in 
terms of blobs in the literature ^ the non-spherical 
character of the blob shape is not explicitly accounted for, 
and it rather is argued that one can characterize the blobs 
by one effective radius ^(r) depending on the radial dis- 
tance r from the cylinder axis. One considers a segmeat 
of the array of length L containing p polymer chains. 'iSl 
On a surface of a cylinder of radius r and length L there 
should then be p blobs, each of cross-sectional area (r) ; 
geometrical factors of order unity are ignored throughout. 
Since the surfaee-area, of the cylindrical segment is Lr, 
we must have liSl 123. ISal 

pe{r) = Lr, CM = (ir/p)i/2 = (r/a)i/2 . (4) 

If the actual non-spherical shape of the blobs (Figure [T]) 
is neglected, the blob volume clearly is of the order of 



FIG. 1: Schematic construction of a blob picture for star polymer (a), and for a bottle brush (b), (c). We assume that along 
the rigid backbone there occur at a regular spacing fa^^ grafting sites where at each grafting site / side chains containing A*' 
effective monomeric units are grafted. If the rigid backbone of the cylindrical bottle brush is oriented along the z-axis, part 
(c) shows a view of the blob partitioning in the xy-plane perpendicular to the backbone, while part (b) shows a cut through 
the cylinder along the xz-plane containing the cylinder axis. A few possible chain configurations are indicated. While for a 
star polymer the blobs have a spherical shape, their radius ^(r) increasing linearly proportional to the distance r from the 
center of the star (a), for the bottle brush the blobs are ellipsoids with three axes ^(r) (in x-direction), proportional to r// (in 
y-direction) , and fo-~^ (in z- direct ion ) , respectively. For a considered chain the x-axis contains the center of mass of the chain. 



^^{r) = (rjuf'l'^. Invoking again the principle that in- 
side a blob self-avoiding walk statistics hold, we have, in 
analogy with Equation ^ , 

i(r) = a[n{T)r , n{r) - K(r)/a]i/- = [r/(aa2)]V2- . 

(5) 

From this result one immediately derives tlie|-«jwep4aw 
decay for the density profile p{r) as follows'^'iSl.l20ll2al 

p{r) = n(r)/e^(r) = a-3[r/((7a2)]-(3i.-i)/2z. 

Using the Flory approximation^ 13 = 3/5 one would 
find p(r) cx r~'^/^ instead. 



Now the average height h of the bottle brush is esti- 
mated by requiring that we obtain all aN monomers per 
unit length in the z-direction along the axis of the bottle 
brush (cf. Figure [T|) when we integrate p{r) from r = 
io r = h, 





This yields, again ignoring factors of order unity, 
7V = (aa)(''-i)/(2.)(;^/^)(.+i)/(2.)^ 



(7) 
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h/a = [aa)^^~-)/i^+-)N^-ni+-) ^ {aaf-^^^ N^-'^ . 

(8b) 

Note that the use of the Flory estimate ly — 3/5 would 
simply yield h cx cr^/^TV'^/^, which happens to be identical 
to the relation that one obtains when one partitions the 
cylinder of height L and radius h into disks of height a^^, 
requiring hence that each chain is confined strictly into 
one such disk. Then each chain would form a quasi-two- 
dimensional self-avoiding walk formed from n'^i^^ = N/n' 
blobs of diameter a~^. Since we have = an"^ again, 
we would conclude that the end-to-end distance R of such 
a quasi-two-dimensional chain is 

= a(aa)-i+3/4"^iV3/4 = aiaa^/^N^^^ . (9) 

Putting then R ~ h, Equation (jSbp results when we use 
there = 3/5. However, this similarity between Equa- 
tions (|8b[) . (ini) is a coincidence: in fact, the assump- 
tion of a quasi-two-dimensional conformation does not 
imply a stretching of the chain in radial direction. In 
fact, when we put the x-axis of our coordinate system 
in the direction of the end-to-end vector R of the chain, 
we would predict that the y-component of the gyration 
radius scales as 

Rgy = a(aa)i/47V3/4 , (iQ) 

since for quasi-two-dimensional chains we have Rg^ oc 
Rgy oc R, all these linear dimensions scale with the same 
power laws. On the other hand, if the Daoud-Cotton- 
like-^1 picture {Figure [T^} holds, in a strict sense, one 
would conclude that Rgy is of the same order as the size 
of the last blob for r = h, 

Rgy = ar = h) = {hlafl^ 

= a(CTa)-'^/(i+'')7V''/(i+'') (11) 

Clearly this prediction is very different from Equa- 
tion (fTU]) . even with the Flory exponent = 3/5 we find 
from Equation that Rgy oc o-^^/^TV^/s^ Hence it is 
clear that an analysis of all three components of the gy- 
ration radius of a polymer chain in a bottle brush is very 
suitable to distinguish between the different versions of 
scaling concepts discussed in the literature. 

However, at this point we return to the geometrical 
construction of the filling of a cylinder with blobs, Fig- 
ure [1] We explore the consequences of the obvious fact 
that the blobs cannot be simple spheres, when we re- 
quire that each blob contains monomers from a single 
chain only, and the available volume is densely packed 
with blobs touching each other. 

As mentioned above, (^(r) was defined from the avail- 
able surface area of a cylinder at radius r, cf. Equa- 
tion Q. It was argued that for each chain in the surface 
area Lr of such a cylinder a surface area ^ (r) is available. 
However, consideration of Figure [1] shows that these sur- 
faces are not circles but rather ellipses, with axes f(J~^ 
and rjj. The blobs hence are not spheres but rather 



ellipsoids, with three different axes: J(J~^ in z-direction 
along the cylinder axes, r// in the y-direction tangential 
on the cylinder surface, normal to both the z-direction 
and the radial direction, and the geometric mean of these 
two lengths (ra~^)^/^, in the radial x-direction. Since the 
physical meaning of a blob is that of a volume region in 
which the excluded volume interaction is not screened, 
this result implies that the screening of excluded volume 
in a bottle brush happens in a very anisotropic way: there 
are three different screening lengths, in the axial z- 

direction, (ra~^^^l'^ in the radial x-direction, and rjj in 
the third, tangential, y-direction. Of course, it remains 
to be shown by a more microscopic theory that such an 
anisotropic screening actually occurs. 

However, the volume of the ellipsoid with these three 
axes still is given by the formula 

(12) 

with ^(r) given by Equation ([5]), and hence the volume 
of the blob has been correctly estimated by the spherical 
approximation. As a consequence, the estimations of the 
density profile pir), Equation ([5]), and resulting brush 
height ft.. Equation (|8b[) . remain unchanged. 

More care is required when we now estimate the linear 
dimensions of the chain in the bottle brush. We now ori- 
ent the X-axis such that the xz-plane contains the center 
of mass of the considered chain. As Figure [T|d indicates, 
we can estimate Rgz assuming a random walk picture in 
terms of blobs. When we go along the chain from the 
grafting site towards the outer boundary of the bottle 
brush, the blobs can make excursions with Az = ±/tT~^, 
independent of r. Hence we conclude, assuming that the 
excursions of the nbiob steps add up randomly 

"blob 

R% = E (/^"')' = "blob(/fT-l)2 , 

2=1 

Rgz = /o-^V'^blob • (13) 

Hence we must estimate the number of blobs nbiob per 
chain in a bottle brush. We must have 

h 

^^blob /" 

"blob = E 1 = y (14) 

Note from Figure [T] that we add an increment 2^(ri) to 
r when we go from the shell i to shell z + 1 in the cylin- 
drical bottle brush. So the discretization of the integral 
in Equation (|14p is equivalent to the sum. From Equa- 
tion (|14p we then find, again ignoring factors of order 
unity 

nbiob = fTi/2fti/2 = (aa)i/(i+'')iV-/(i+'^) , (15) 
and consequently 
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With Flory exponents we hence conclude ribiob oc 
^5/8^3/8^ and thus Rg^ oc /cr-^/i^iV^/ie. The estima- 
tion of Rgy is most delicate, of course, because when we 
consider random excursions away from the radial direc- 
tions, the excursions Ay — iC(^) S'^'^ non-uniform. So 
we have instead of Equation (I13|) . 

h 

^^blob /. 

1=1 

in analogy with Equation This yields Rgy — 

^-1/2^3/2 g^j-^i^ hence 

= a-l/4/i3/4 = ^(^„)-(2.-l)/(2.+2)^3./(2.+2) 

while the corresponding result with Flory exponents 
is Rgy CX cr^i/i^iV^/is. These different results for 
Rgx OC /i {Equation ((8b)) }. Rgy {Equation (flSjl } and 
-Rgz {Equation (|16p } clearly reflect the anisotropic struc- 
ture of a chain in a bottle brush. Note that the result 
for Rgy according to Equation (fT5|) is much larger than 
the simple Daoud-Cotton-like prediction. Equation pTjl . 
but is clearly smaller than the result of the quasi-two- 
dimensional picture. Equation pn]) . 

As a consistency check of our treatment, we note that 
also Equation ([7]) can be formulated as a discrete sum 
over blobs 

h h 


(19) 

which yields Equation ([5a|) . as it should be. 

Finally we discuss the crossover towards mushroom be- 
havior. Physically, this must occur when the distance 
between grafting points along the axis, fa~^, becomes 
equal to aN" . Thus we can write 

aN^hiaaN""), (20) 

where we have absorbed the extra factor / in the scaling 
function h{C,). We note that Equation (|8b| results from 
Equation (I^Dl) when we request that 

/i(C> 1) oc C^i-'')/(i+''), (21) 

and hence a smooth crossover between mushroom be- 
havior and radially stretched polymer conformations oc- 
curs for aaN" of order unity, as it should be. Analogous 
crossover relations can be written for the other linear di- 
mensions, too: 

Rg. = aN'^Rg.iO, RgAC » 1) (X ^-(2.+i)/(2.+2)^ 

(22) 

Rgy = aN-^RgyiO, Rgy{C » l) CX ^ " ^ D / (2^ + 2) , 

(23) 



The agreement between Equation ([T6|) and Equation ([22]) 
or Equation and Equation (^5]) . respectively, pro- 
vides a check on the self-consistency of our scaling argu- 
ments. 

However, it is important to bear in mind that 
the blob pieture of polymer brushes, as deveioped by 
Alexander,!^ Daoud and Cotton,*! Halperin'^'^ and 
many others, is a severe simplification of reality, since 
its basic assumptions that (i) all chains in a polymer 
brush are stretched in the same way, and (ii) the chain 
ends are all at a distance h from the grafting sites, 
are not true,— ■'fceiatrparita |ba,sed|-ejii— the self-consistent 
field theor^ M imUpJaJp^ ^ .-^nd computer 
simulationstal2dMl2cll72l,yMy,yMy ^ave shown 

that chain ends are not confined at the outer bound- 
ary of the brush, and the monomer density distribution 
is a nontrivial function, that cannot be described by the 
blob model. However, it is widely believed that the blob 
model yields correctly the power laws of chain linear di- 
mensions in terms of grafting density and chain length, so 
the shortcomings mentioned above affect the pre-factors 
in these power laws only. In this spirit, we have extended 
the blob model for brushes in cylindrical geometry in the 
present section, taking the anisotropy in the shape of 
the blobs into account to predict the scaling behavior 
of both the brush height h (which corresponds to the 
chain end-to-end distance R and the x-component Rgx 
of the gyration radius) and of the transverse gyration ra- 
dius components Rgy, Rgz- To our knowledge, the latter 
have not been considered in the previous literature. 



III. Monte Carlo Methodology 

We consider here the simplest lattice model of poly- 
mers under good solvent conditions, namely, the self- 
avoiding walk on the simple cubic lattice. The backbone 
of the bottle brush is just taken to be the z-axis of the 
lattice, and in order to avoid any effects due to the ends 
of the backbone, we choose periodic boundary conditions 
in the z-direction. The length of the backbone is taken to 
be Liy lattice spacings, and the lattice spacing henceforth 
is taken as one unit of length. Note that in order to avoid 
finite size effects Lf, has to be chosen large enough so that 
no side chain can interact with its own "images" created 
by the periodic boundary condition, i.e. Rgz <C L},. 

We create configurations of the bottle brush poly- 
mers appij4ii& jthfi-^runed-enriched-Rosenbluth method 
(PERM).l^l^'l^l42, 87, 88, 89, 90] rpj^jg algorithm is 

based on the biased cbain growth algorithm of Roscn- 
bluth and Rosenbluth,'^! and exteads it by re-sampling 
technique ("population control"'^) with depth-first 
implementatjiep. Similar to a recent study of star 
polymers'^ 1^ all side chains of the bottle brush are 
grown simultaneously, adding one monomer to each side 
chain one after the other before the growth process of 
the first chain continues by the next step. 
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FIG. 2: Snapshot picture of a bottle brush polymer with 
1/6 — 128 (7 = 1/4, and A'' = 2000 on the simple cubic lattice. 
Note that diflerent chains are displayed in different colors to 
distinguish them, and the periodic boundary condition is un- 
done for the sake of a better visibility in the visualization of 
this configuration. 



When a monomer is added to a chain of length n — 1 
(containing n monomers) at the nth step, one scans the 
neighborhood of the chain end to identify the free neigh- 
boring sites of the chain end, to which a monomer could 
be added. Out of these rifrec sites available for the addi- 
tion of a monomer one site is chosen with the probability 
Pn^i for the zth direction. One has the freedom to sam- 
ple these steps from a wide range of possible distribu- 
tions, e.g. pn^i — 1/nfroc, if one site is chosen at random, 
and this additional bias is taken into account by suitable 
weight factors. The total weight Wn of a chain of length 
n with an unbiased sampling is determined recursively 
by Wn = njJ^^^Wfc = Wn-iWn- While the weight Wn is 
gained at the nth step with a probability p„^i, one has to 
use Wn/pn,i instead of Wn- The partition sum of a chain 
of length n (at the nth step) is approximated as 

M„ 

Z„ « Z„ = M- 1 ^ W„ (a) , (24) 

a=l 

where M„ is the total number of configurations {a}, and 
averages of any chain property (e.g. its end-to-end dis- 
tance, gyration radius components, etc.) A(a) are ob- 
tained as 

M„ 

J _ ^ a=l /ocN 



As is weU-known from Ref. IM], for this original biased 
sampling's^ f}iQ statistical errors for large n are very 
hard to contsal. This problem is alleviated by popula- 
tion control. 'S3 One introduces two thresholds 

W+^C+Z„, W~^C-Zn , (26) 



where Zn is the current estimate of the partition sum, 
and C_|_ and C_ are constants of order unity. The opti- 
mal ratio between C+ and C_ is found to be C+ / C_ ^ 10 
in general. If Wn exceeds W^ for the configuration a, 
one produces k identical copies of this configuration, re- 
places their weight Wn by Wn/k, and uses these copies 
as members of the population for the next step, where 
one adds a monomer to go from chain length n to n -I- 1. 
However, if Wn < W~ , one draws a random number rj, 
uniformly distributed between zero and one. If 77 < 1/2, 
the configuration is eliminated from the population when 
one grows the chains from length n to n -I- 1 . This "prun- 
ing" or "enriching" step has given the PERM algorithm 
its name. On the other hand, if 77 > 1/2, one keeps the 
configuration and replaces its weight Wn by 2Wn- In a 
depth-first implementation, at each time one deals with 
only a single configuration until a chain has been grown 
either to the end of reaching the maximum length or to 
be killed in between, and handles the copies by recursion. 
Since only a single configuration has to be remembered 
during the run, it requires much less memory. 

In our implementation, we used Wn = 00 and W^ = 
for the first configuration hitting chain length n. For the 
following configurations, we used Wn — CZn(cn/ cq) and 
Wn — 0.15W+, here C = 3.0 is a positive constant, and 
c„ is the total number of configurations of length n al- 
ready created during the run. The bias of growing side 
chains was used by giving higher probabilities in the di- 
rection where there are more free next neighbor sites and 
in the outward directions perpendicular to the backbone, 
where the second part of bias decreases with the length 
of side chains and increases with the grafting density. To- 
tally 10^ ^ 10^ independent configurations were obtained 
in most cases we simulated. 

Typical simulations used backbone lengths L;, = 
32, 64, and 128, / = 1 (one chain per grafting site of the 
backbone, although occasionally also f — 2 and / = 4 
were used), and grafting densities a = 1/32, 1/16, 1/8, 
1/4, 1/2 and 1. The side chain length TV was varied up to 
N = 2000. So a typical bottle brush with Lb ^ 128, a = 
1/4 (i.e., nc = 32 side chains) and N — 2000 contains 
a total number of monomers A^tot — Li, + ncN = 64128. 
Figure [2] shows a snapshot configuration of such a bottle 
brush polymer. Note, tiiajt-paest other simulation algo- 
rithms for polymers'S3''S3.lSa|86[ -^vould fail to produce a 
large sample of well-equilibrated configurations of bottle 
brush polymers of this size: dynamic Monte Carlo al- 
gorithms using local moves involve a relaxation time (in 
units of Monte Carlo steps per monomer [MCS]) of or- 
der N'^ where z = 2j/ -I- 1 if one assumes that the side 
chains relax independieHtly of each other and one apr 
plies the Rouse model '23' in the good solvent regime. 'Sa 
For N — 2000 such an estimate would imply that the 
time between subsequent statistically independent con- 
figurations is of the order of 10^ MCS, which clearly is 
impractical. While the pivot algorithm'^ would pro- 
vide a significantly faster relaxation in the mushroom 
regime, the acceptance rate of the pivot moves quickly 



deteriorates when the monomer density increases. This 
algorithm could equilibrate the outer region of the bottle 
brush rather efficiently, but would fail to equilibrate the 
chain configurations «ea*,the backbone. The configura- 
tional bias algorithmic I^J would be an interesting alter- 
native when the monomer density is high, but it is not 
expected to work for very long chains, such as = 2000. 
Also, while the simple enpchment technique is useful to 
study both-Sitar polymers'2^ and polymer brushes on flat 
substrates'SSl it also works only for chain lengths up to 
about N — 100. Thus, existing Monte Carlo simulations 
of one-component bottle brushes under good solve^l con- 
ditions either used the bond fluctuation model- on 
the simple cubic lattice appk'ing local raoves with side 
chain lengths up to iV = 4l[23 or TV = 64129 or the piyet 
algorithm with side chain lengths up to = 80, 
but considering flexible backbone of length L = SQO. An 
alternative approach was followed by Yethiraj'^y who 
studied an off-lattice tangent hard sphere model by a 
pivot algorithm in the continuum, for N < 50. All these 
studies addressed the question of the overall conforma- 
tion of the bottle brush for a flexible backbone, and did 
not address in detail the conformations of the side chains. 
Only the total mean square radius of gyration 
side chains was-estimated occasicpally, obtainingl^i 
Rl (X iVi-2 oM i?2 cx andlai cx N^-^. How- 

ever, due to the smallness of the side chain lengths used 
in these studies, as quoted above, these results have to 
be considered as somewhat preliminary, and also a sys- 
tematic study of the dependence on both N and a was 
not presented. We also note that the conclusions of the 
quoted papers are somewhat contradictory. 

An interesting alternative simulation method to the 
Monte Cada study of polymeric systems is Molecular 
Dynamics, of course. While in corresponding stud- 
ies of a— bead-spring model of polymer chains for flat 
brushes chain lengths N up to N — 200 were used, 
for chains grafted to thin cylinders the three chain 
lengths N = 50, 100 and 150 were used. For TV-p 50, 
also four values-,of grafting density were studied, Mu- 
ral and Grest'23 used these data to test the scaling pre- 
diction for the density profile, Equation but found 
that p{r) is better compatible with p{r) oc r~°-^ rather 
than p{r) cx r""-^^. However, for = 50 the range where 
the power law is supposed to apply is very restricted, and 
hence this discrepapey was not considered to be a prob- 
lem for the theory. 123 

Thus, we conclude that only due to the use of the 
PERM algorithm has it become possible to study such 
large bottle brush polymers as depicted in Figure [H Nev- 
ertheless, as we shall see in the next section, even for such 
large side chain lengths one cannot yet reach the asymp- 
totic region where the power laws derived in the previous 
section are strictly valid. 



IV. Monte Carlo Results for One- 
Component Bottle-Brush Polymers 
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FIG. 3: Log-log plot of the mean square end-to-end distance 
components {R^) (a), {Ry) (b), and {R^) (c) versus side chain 
length A'^, for various choices of the grafting density a as indi- 
cated. All data refer to / = 1 (one chain per possible grafting 
site) and TV > 5. Note that the x-direction for every chain is 
the normal direction from its center of mass position to the 
bottle brush backbone. All data are for Lb = 128. The chain 
mean square linear dimensions are all normalized by N^'^ . 



In Figures [3] and [H we present our data for the mean 
square end-to-end radii and gyration radii components of 
the side chains. Note that for each chain configuration a 
separate local coordinate system for the analysis of the 
chain configuration needs to be used; while the z-axis is 
always oriented simply along the backbone, the x-axis 
is oriented perpendicular to the z-axis and goes through 
the center of mass of the chain in this particular config- 
uration. The y-axis then also is fixed simply from the 
requirement that it is perpendicular to both the x- and 
z-axes. 

For a polymer mushroom (which is obtained if the 
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FIG. 5: Log-log plot of Ri ^ Rl + divided by iV^" 
vs. ( — aN" , including all data for N > 5, and three choices of 
Lt as indicated (a), or alternatively removing data affected by 
the finite size of the backbone length via the periodic bound- 
ary condition (b). The slope indicated in (b) by the dotted 
straight line corresponds to the scaling estimate from Equa- 
tion 2(1 - v)/{l + v)^ 0.519. 
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FIG. 4: Log-log plot of the mean square gyration radius com- 
ponents (Rgx) (a), {Rgy) (b), and {R?gz) (c), versus side chain 
length TV. Only data for A*' > 10 are included. All data 
are for / = l,Lb — 128, and various choices of a. All chain 
mean square linear dimensions are normalized by N^'^ with 
i/ = 1/3 « 0.588. 



grafting density a is sufficiently small) we expect that all 
chain linear dimensions scale as iV", for sufficiently long 
chains. Therefore we have normalized all mean square 
linear dimensions in Figuses H,and S] by a factor N'^'^, 
using the theoretical value'^'l^ v = 0.588. However, as 
we see from Figures |31 H] in the range 10 < iV < lO'^ dis- 
played there, even for the smallest a presented, where 
a single side chain is simulated, the shown ratios are 
not coHStaat. This fact indicates that corrections to 
scalinglS^l^ should not be disregarded, and this prob- 
lem clearly complicates the test of the scaling predic- 
tions derived above. For the largest value of a included 
(cr = 1/4), the irregular behavior of some of the data 
(Figure ^jp,c, Figure |3)d) indicate a deterioration of sta- 



tistical accuracy. This problem gets worse for increasing 
number of side chains ric- 

A further complication is due to the residual finite size 
effect. Figure [5] shows a plot of R\ = + Ry vs. the 
scaling variable ( = o-N". One can recognize that for 
small Lb but large N and not too large cr systematic 
deviations from scaling occur (Figure [5^) , which simply 
arise from the fact that an additional scaling variable 
(related to (i?^)^/^/Lf,) conies into play when (i?^)^/^ no 
longer is negligibly smaller than Li,. While for real bot- 
tle brush polymers with stiff backbone effects due to the 
finite lengths of the backbone are physically meaningful 
and hence of interest, the situation is different in our 
simulation due to the use of periodic boundary condi- 
tions. The choice of periodic boundary conditions is mo- 
tivated by the desire to study the characteristic structure 
in a bottle brush polymer far away from the backbone 
ends, not affected by any finite size effects. However, if 
(Rg^)^^^ becomes comparable to Lt, each chain interacts 
with its own periodic images, and this is an unphysical, 
undesirable, finite size effect. Therefore in Figure [SJd, the 
data affected by such finite size effects are not included. 
One finds a reasonable data collapse of the scaled mean 
square end-to-end distance when one plots the data ver- 
sus the scaling variable C = o-N^. These results consti- 
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FIG. 6: Log-log plot of {R%)/N'^'' (a), {Rly)/N'^'' (b) and 
{R^gz) /N"^" (c) versus the scaling variable C, = aN'^ , using 
the data in Figure |4] but for A*' > 50. The slopes indicated 
by dashed straight lines illustrate the scaling estimates from 
Equation ((21]), 2(1 - -f !/) « 0.519 (a), from Equa- 

tion (1231), -{2u - l)/(v + 1) ?a -0.111 (b), and from Equa- 
tion (EU, -(21/ + l))/(y +1)^ -1.370 (c). 



tute the first comprehensive test of the scalmg relations 
for bottle brush polymers, Equations ((2T|) - ((23l) , and the 
consistency between the data and the proposed scaling 
structure in terms of the variable C is indeed very grati- 
fying. On the other hand, it is also evident from FigureslH] 
and [6] that only a mild stretching of the side chains away 
from the backbone is observed, and in this region one is 
still far away from the region of strong stretching, where 
the simple power laws Equations ((8b)) . ((T5)) and ((TB)) ap- 
ply. Obviously, the crossover from the simple mushroom 
behavior (observed for C = uA^'^ <C 1) to the strongly 
stretched bottle brush (observed for C ^ 1) takes at least 
one decade of the scaling variable C,. There is a rather 
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FIG. 7: Log-log plot of the mean-square gyration radius com- 
ponents {R%,,n,) / N^^"^ (a), {Rly ,,,^) / N^^"^ (b) and {R%^,,,,) 
(c) versus side chain length A'^ with U2 = 3/4. As in Fig- 
ure (4] only data for N > 5 are included. Various choices of 
a are included. Note that the walls of the slit are located at 
z — ±((T~^ -|- l)/2 and are strictly repulsive hard walls. One 
chain end is fixed at the x-axis at z = 0, and since the x-axis 
represents the backbone if this polymer in a slit is taken as a 
model of a disk-shaped section of a bottle brush, all sites of 
the X-axis are excluded from occupation of the monomers of 
the (side) chain as well. 



gradual and smooth crossover rather than a kink-like be- 
havior of the scaling function. While in the x-component 
at least near ^ = 10 a weak onset of stretching can be 
recognized, hardly any evidence for the contraction of the 
y and z-components is seen. 

In order to test what one would expect if the picture 
of the chains as quasi-two-dimensional self-avoiding walks 
were correct, we have also studied single chains grafted to 
a straight backbone of length a, a = 1, 1/3, 1/5, . . . 1/63 
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FIG. 8: Log-log plot of {Rg^,siit) / N^'' (a), {R^ 
^gzEht) /^^'^ (c) versus the scaling variable 



Ht)/iV^^ and 
= criV. Only 

data for A*' > 10 are included. The slopes indicated by dashed 
straight lines illustrate the scaling exponent implied by Equa- 
tion (|27[) . namely 3/2u — 2 (a,b) and —2 (c), respectively. In 
(c) is replaced by + 1 to remove a finite-size effect. 



which are confined between two parallel repulsive infinite 
walls. The grafting site of the chains was chosen at the 
site located symmetrically between the confining walls in 
the slit {if the grafting site is chosen to be the origin of 
the coordinate system, the confining hard walls occur at 
z = ±(tT~^ -I- l)/2}. Figure [7] gives log-log plots of the 
gyration radii components of such chains confined to such 
disk-like slits, which we denote as Rga, slit, a = x,y and 
z, in order to distinguish then from the actual gyration 
radii components of the side chains in a bottle brush 
polymer. One can clearly see that (i?g^,siit), (^gy.siit) 
scale as N'^^'^ = N'^/'^ for large N, as it must be. These 
data are very similar to data for chains confined between 
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FIG. 9: Log-log plot of {R%) KR%,,ii,) (a), {R. 
(b) and {R%) / {Rlz.sut) versus A'', for values of a strictly cor- 
responding to each other, as indicated. 



repulsive "W^Us without a grafting to a piece of a rigid 
backbone, '23 of course. In the corresponding scaling plot 
(Figure [5]) one can see that both 7?^^ ^ and have 
the simple crossover scaling behavior^. Ii~ 



gy, siit 



(27) 

with /^(C) cx fyiC) cx (2(1.2-1')/!^ ^ith 1^2 =rA/A as ex- 
pected, and seen in related previous work. 12211221 Obvi- 
ously, although we use in Figures [7] and [5] the same range 
of N and a as in Figures [J] and [51 the behavior is rather 
different. As expected from our scaling analysis pre- 
sented above, the confinement that a chain experiences 
in a bottle brush due to the presence of the other chains 
is much weaker than the confinement of a chain in the 
equivalent disk-like sector between confining walls. This 




FIG. 10; Radial distribution function p{r) plotted vs. r, for 
N = 500 (a), TV = 1000 (b), and N = 1500 (c), for various 
choices of cr, as indicated. Note that due to the discreteness 
of the lattice, the number A'^(r) of monomers in the interval 
[r, r + dr] is not normalized by the factor nr that applies in the 
continuum limit, but by the number Nr of lattice sites (x,y) 
satisfying the constraint = +y'^, i.e. N(r) = N and 
p{r) = N{r)/Nr. All data refer to U = 128 and / = 1. 



fact is demonstrated very directly in Figure [51 where the 
ratios 

,siit) plotted vs. N for the corre- 
sponding values of a. If the hypothesis of quasi-two di- 
mensional confinement were valid, we would expect these 
ratios to be constants. Obviously, this is not the case. 

As a final point of this section, we discuss the distri- 
bution of monomers (Figures (TU] and llip and chain ends 
(Figures [T^ and [T5)) in the simulated model for the bottle 
brush polymer. Unlike corresponding radial density dis- 
tributions for the off-lattice bead-spring model of Murat 
and Grest,l23 where for small distances close to the back- 
bone a kind of "layering" was observed, we see a rather 



FIG. 11: Density distribution function p{z) plotted vs. the 
coordinate z along the backbones, for A'' — 500 (a). A'' = 1000 
(b), and N = 1500 (c), and various choices of a, as indicated. 
This distribution is normalized by choosing p{z) = 1 



smooth behavior also for small distances (Figure [10]). For 
larger distances, the behavior is qualitatively very sim- 
ilar. Again we fail to provide a clear-cut evidence for 
the predicted power law decay. Equation Note, how- 
ever, that this power law is supposed to hold only in 
the strong stretching limit, where Equation (|8bp is ob- 
servable (which we do not verify either), and in addition 
the stringent condition 1 <C r ^ /i needs to be obeyed, 
to have this power law. Although our simulations en- 
compass much longer chains than every previous work 
on bottle brushes, we clearly fail to satisfy this double 
inequality. 

Turning to the distribution along the backbone (Fig- 
ure [11]), the periodicity due to the strictly periodic spac- 
ing of grafting sites is clearly visible. If desired, one could 
also study a random distribution of grafting sites along 
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the backbone, of course, but we have restricted attention 
to the simplest case of a regular arrangement of grafting 
sites only. The distribution of chain ends p{z) exhibits 
an analogous periodicity for small values of cr only, while 
Peiz) is approximately constant (Figure I12p for larger 
values of a. 

Most interesting is the radial distribution of chain ends 
fFigure fO)) . One can see an increasing depression of Pe (r) 
for small r when a increases. Again these dateiare similar 
to the off-lattice results of Murat and Grest.'231 In no case 
do we see the "dead zone" predicted by the self-consistent 
field theory in the strong stretching limit, "^^'-^^ however. 



V. Phase Separation in Two-Component 
Bottle Brushes: Theoretical Background 



We still consider a bottle brush polymer with a strictly 
rigid straight backbone, where at regularly distributed 
grafting sites (grafting density a) side chains of length 
N are attached, but now we assume that there are two 
chemically different monomeric species, A and B, com- 
posing these side chains with = Nb — N. These 
systems have found recent attention, suggestHier-tii|e-*os- 
sibility of intramolecular phase separation. l5a|53.|58|l jj-^ 
a binary system pairwise interaction energies e^A, ^ab, 
and esB are expected, and consequently phase separation 
between A and B could be driven by the Flory-Huggins 
parameter 

X = Zce/kBT, e = EAB - i^AA + ^bb)/^ (28) 

where Zc is the coordination number of the lattice, like in 
the Flory-Huggins lattice model of phase separation in 
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FIG. 14: Schematic description of perfect phase separation of 
side chains in a cylindrical binary (A,B) bottle brush poly- 
mer into a "Janus cylinder" structure, where the A-chains 
occupy the lower hemicylinder and the B-chains occupy the 
upper hemicylinder. The interface between B-rich and A-rich 
phases is assumed to be oriented in the yz-plane. The figure 
indicates the characterization of this order in terms of a local 
order parameter, namely the vector 'il){z) oriented normal to 
the interface at every grafting site. The absolute value of this 
vector characterizes a suitable average of the local phase sep- 
aration in a disk of the cylinder located at z (see text). For 
perfect order {| '(/'(•^) 1) = 1 ^-nd tp{z) is oriented along the 
same axis for the whole bottle brush polymer. 



a binary polymer blend. B S EE lil E3 However, since 
the side chains are grafted to the backbone, only in- 
tramolecular phase separation is possible here. Actu- 
ally, the energy parameter e suffices for dense polymer 
blends or deHS«-|bkj(sk,oe|)olymer melts, where no solvent 
is present, I4a ISfl ISlI, l5j Pqj. ^jjg problem of intramolec- 
ular phase separation in a bottle brush, the solvent can- 
not be disregarded |-an|d-,hence the enthalpy of mixing 
rather is written as.l^l^ 



— [ XASfj^AirJcj^sir) + XBS<l>Bir)(t)s{r) 
+ XAB0A(r)0B(r)], (29) 

where v is the volume per monomer, 0^(r^, (j)B{r) and 
(j)s ('^) are the local volume fractions of monomers of types 
A and B, and the solvent density, respectively, and three 
pairwise interaction parameters Xas,XbSi and xab en- 
ter. The latter xab corresponds to the x-parameter writ- 
ten in Equation while the former two control the 
solvent quality for polymers A and B, respectively. Of 
course, using the incompressibility condition 



(r) + (t)B [r) + (t>s [r) = 1 



(30) 



the solvent density can be eliminated from the problem 
(but one should keep in mind that in the free energy 
density there is an entropy of mixing term (ps (r) In <j)s {r) 
present ISal). In the framework of the lattice model stud- 
ied here, v = 1 and solvent molecules are only implicitly 
considered, identifying them with vacant sites. For sim- 
plicity, the following discussion considers only the most 
symmetric case, where xas = Xbs, and the numbers of 
A chains and B chaips are equal. 

Stepanyan et al. used Equations P5|) and (PO)) . 
as the starting point of a self-consistent field calcula- 



tion, adding conformational free energy contributions ac- 
counting for the entropy associated with the stretching 
of Gaussian chains. It is found that when xab exceeds 
a critical value X*AB^ intramolecular phase separation of 
"Janus cylinder" -type occurs (Figure [T4|) . Then an in- 
terface is formed, containing the backbone of the bottle 
brush, such that there is an excess of A-monomers below 
the interface and an excess of B-mononiers above the in- 
terface (or vice versa). Assuming that at the position z of 
the backbone the interface is oriented in x-direction, we 
can describe this Janus-type phase separation in terms 
of an one-dimensional order parameter density 

V'(z) = {n+-n\ + n-]^-n-g)/{n\ + n-^ + 7i++ng), (31) 

where and are the numbers of A(B) monomers 
in the interval -I- dz] with x > 0(nj^,n^) and x < 
0(n^,njj), respectively. Since we shall see below that 
the orientation of the interface is an important degree of 
freedom, we may consider V'(^) as the absolute value of 
a vector order parameter ipiz), and the direction of ip{z) 
is defined such that takes a maximum. However, 

Stepanyan et al. did not consider the possibility of 
an inhomogeneity along the z-axis, and they also did not 
derive how the order parameter depends on the parame- 
ters N, (7, and the various x paranieters of the problem 
{Equation ((29l) }. Stepanyan et al.l^ also assume that 
the distribution of chain ends is a delta function at the 
radius (or "height" h) respectively) of the bottle brush, 
where h is given by /i/a = {cray/^N^/^ in the good sol- 
vent regime. They then find that phase separation occurs 
for 



(32) 



X*AB « yVN, 



but argue that this result holds only for a "marginal sol- 
vent" rather than a good solvent, due to the mean- field 
approximation used which neglects that inside a blob 
all binary contacts are avoided, estimating the number 
of contacts simply proportional to the product (f>A<t'Bi 
cf. Equation (|29p . and neglecting the correlations due 
to excluded volume. The regime of marginal solvents is 
reached near the 0-point (which o©|;urs for 1 — 2xas — 
1 — 2xbs — 0), and requires that ^ 



< 1 - 2xAS < N-^/^ 



(33) 



where again prefactoss of order unity are omitted. 

Stepanyan et al. extend Equation ([5^ by a simple 
scaling- type argument, stating that the transition from 
the mixed state to the separated state occurs when the 
energy of the A-B contacts per side chain, AE, is of the 
order of fc^T. They estimate this energy as 



AE/keT^ Np{^)XAB, 



(34) 



where p(0) is the probability of the A-B contact, 
which depends on the average volume fraction (j) of the 
monomers inside the brush. According to the mean field 
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FIG. 15: Same as Figure 1141 but allowing for a long wave- 
length variation of the vector ^p{z) characterizing the local 
interface orientation in the Janus cylinder. 
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FIG. 16: Schematic sketch of the temperature dependence of 
the correlation length ^ of the ferromagnetic XY-model in a 
cylinder of radius R. 



theory, oc 0, where (f) — N/h'^, h being the ra- 

dius of the cyhndrical brush. Using h cx N^^^ (Equa- 
tion (jSbp ) one finds p{(j)) (x N'^^^ and using this result 
in Equation ((34|) implies AE/ksT cx N^^^xab, and from 
AE/ksT = 1 we recover Equation ([32]) . 

The merit of this simple argument is that it is readily 
extended to other cases: e.g., for a 6'-solvent,|-we have^i^ 



h (X 7V^/^ and hence p{(f>) oc N yieldi 



mg 



m 



X*AB « N-^^^ , e - solvent. 



(35) 



For the poor solvent case, the bottle brush should col- 
lapse to a cylinder densely filled with monomers, and 
hence h oc N^^"^, ^ = 1, and thus 



Xab N ^ , poor solvent. 



(36) 



Note that Equation is the same relation as for 

a dense bult-DeipaBP-, blend or block copolymer melt, 
respectively. I^l^l5iy,y 

The situation is more subtle in the good solvent case, 
however, since there the probability of contact is no 
longer given by the mean field result p{4>) c>c (j) but rather 

by ad 



p{4>) oc ^5/4 cx , 



(37) 



where the Flory approximation = 3/5 was again used, 
with (j) oc iV"i/2 when h oc N^/*. Equations and 



37)1 now imply 



Xab ^ "^^^ : good solvent 



(38) 



However, these crude estimates do not suffice to prove 
that for XAB > Xab ^^'^^ lo'^S range order of this "Janus 
cylinder" type (Figure fT4|) is established. Hsu et al.'^ 
suggested that there is also the need to consider varia- 
tions of the direction of the order parameter ipiz) along 
the z-direction (Figure [15]). It was argued that for any 
finite side chain length N also the cylinder radius (or 
brush "height") h is finite, and hence the system is 
one-dimensional. In one-dimensional systems with short 
range interactions-gtt peazero temperatures no long range 
order is possible. 'i^SlliSiJ xhe situation depicted in Fig- 
ure [15] is reminiscent of the one-dimensional XY-model 
of a chain of spins on a one-dimensional lattice where 
each spin at site i is described by an angle (pi in the xy- 
plane, with < i < 27r, and where neighboring spins are 
coupled. This coupling is described by the Hamiltonian 



H = — J^C0s((/3j 



(39) 



when Si — (cos ipi, sin tpi) is a, unit vector in the xy-plane. 
While mean field theory predicts that ferromagnetic or- 
der occurs along the chain, for ferromagnetic exchange 
J > and temperatures T less than the critical tem- 
perature T*^^ which is of the order of J/A:p7-|tiie-,exact 
solution of this model. Equation showstiSa lifli that 
Tc = 0, since ferromagnetic long range order is unstable 
against long wavelength fluctuations. One can show that 
the spin-spin correlation function for large z = a{j — i) 
decays to zero exponentially fast. 



{Si ■ Sj) cx exp[-z/^], z 



(40) 



where a is the lattice spacing of this spin chain. The 
correlation length ^ gradually grows as the temperature 
is lowered, 



^ = 2a( J/fcsT) 



(41) 



Equation (|4ip is at variance with mean field theory, which 
rather would predict '^2211 (the index MF stands for "mean 
field" throughout) 



MF 



1) 



I^MF = 1/2 



(42) 



This critical divergence at a nonzero critical tempera- 
ture r*^^ is completely washed out by the fluctuations: 
rather than diverging at T*^^, the actual correlation 
length ^ {Equation (jH])} at T^'^ still is only of the order 
of the lattice spacing. 

This consideration can be generalized to spin sys- 
tems on lattices which have a large but finite size in 
(d — 1) dirfteiis4©HS and are infinite in one space dimen- 
sions only. 'i23> E.g., when we consider an infinitely 
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long cylinder of cross section irR'^ we expect that Equa- 
tions (HpTi ' ^ I are replaced by the finite size scaling 
relationllQira 



e = ?oo |(i?/eoo) , |(C»1) = 1, |(C«1)«C, (43) 

where ^oo is the correlation length of the XY model on a 
lattice which is infinitely large in all d = 3 directions of 
space. 



Coo(x(T/re-l)- 



(44) 



^ = 2TTT{T)R^/kBT w 2tt [R^/a){J/kBT), T ■ 



where again ksTc oc J (but with a smaller constant of 
proportionality than that in the relation ksT^^^ oc J), 
and I' is the corpeiation length exponent of the XY 
model {ly « 0.67).'^^ Equation (|43|) describes a smooth 
crossover of the ferromagnetic correlation length describ- 
ing spin correlations along the axis of the cylinder from 
bulk, Equation (|44p . to a quasi-one-dimensional varia- 
tion. For T ^ X—tlae-correlation length resembles Equa- 
tion Ell), sinceliOa'lifld 

0. 

(45) 

the "helicity modulus" (also called "spin wave stiffness") 
r(T) is of order J for T ^ and shows a critical vanish- 
ing as T ^ Tj, from below, in the thermodynamic limit 
i? — !■ oo. However, for finite R there is a finite size round- 
ing of this singularity of r(r TJ, such that r(r = 
is of order l/R, and hence a smooth crossover to Equa- 
tion (|43| occurs near Tc. Figure [16] gives a qualitative 
account of this behavior. For more detaiis-|Of-yiis finite 
size crossover we refer to the literature. 'i^SliSJ g^t we 
suggest a qualitatively similar behavior for the domain 
size ^ of segregated A-rich and B-rich domains in binary 
bottle brush polymers. So, when we study the correla- 
tion function of the order parameter considered in Equa- 
tion dSll) 



G^{z) = {il^{z') ■ 'ijj{z' + z)) oc exp(-z/^^) 



(46) 



we expect that the correlation length describing the lo- 
cal phase separation in the direction along the backbone 
of the bottle brush polymer remains of order unity as 
long as Xab exceeds x^^^ distinctly. For xab near Xab 
we expect that becomes of order h, the radius of the 
bottle brush. For xab ^ X*abi expect oc h?XAB, 
by analogy with Equation ([15)) . Unfortunately, the test 
of those predictions even in the poor solvent case where 
Xab^ is largest, is rather difficult, since the prefactor in 
the relation x*ab ^ 1/-^ is unknown. 



VI. Monte Carlo Results for Binary Bottle 
Brush Polymers 

We use the same lattice model as considered before 
in our Monte Carlo study of the chain conformations of 



one-component bottle brush polymers, but now in the 
construction of the weights Wn{cy) we have to take into 
account that the partition function now is 



^ - 2^1 Iab 



(47) 



with (remember that we restrict attention to the choice 

f^AA ~ £Bb) 

q = cxp{-eAA/kBT), qab = ex.p{~eAB /kBT) . (48) 

In Equation (|47p the numbers of non-bonded occupied 
nearest-neighbor monomer pairs AA, BB and AB are de- 
noted as rnAAjWBB and mAB, respectively. Note that 
the sum in Equation ([17[) extends over all possible con- 
figurations {a} of the bottle brush polymer. The choice 
Q = QAB = 1 corresponds to the previously studied one- 
component bottle brush under good solvent conditions, 
while the choice q — qab > 1 corresponds to variable 
solvent quality for the one-component brush (note that 
q — QAB means x = 0, Equation ([28]) . and also xab 
which is proportional to x then vanishes: this means 
there is no chemical incompatibility between A and B 
any longer, no physical difference between A and B pr- 
ists any more). From previous work on single chains'^ 
we know that the 0-point occurs for qg = 1.3087. There- 
fore we varied q in the range 1 < q < 1.5; q — 1.5 hence 
falls in the regime of poor solvent quality already. Of 
course, in order to have rather compact configurations of 
cylindrical bottle brushes a choice of much larger q would 
be desirable. However, the efficiency of the PERM algo- 
rithm quickly deteriorates with increasing q: for q — 1.5 
we encounter already for rather small values of the side 
chain length N such as = 18 and a backbone length 
of Lb = 64 huge statistical fluctuations. The total size 
of the bottle brush polymer under poor solvent condi- 
tions reached here, A^tot — NaLi, + Li, = 1216, is al- 
most two orders of magnitude smaller than the maxi- 
mal size studied under good solvent conditions! How- 
ever, all known simulation algorithms for polymers suffer 
from difficultieSj-ef leauilibration in the limit of very dense 
configurations. lSJ.lfia,lSal 

Since we are mostly interested in the high grafting 
limit ((T = 1, so the number of side chains ric = 
aLt = Lb) in the PERM algorithm where all side chains 
grow simultaneously we use a bias factor such that side 
chains are grown with higher probability in the direc- 
tions perpendicular to the backbone. This additioBaJ 
bias (which is not present in the standard Rosenbluth 121 
and PERM-'*-l methods) must be taken into account by 
suitable weight factors. About 10^ independent configu- 
rations were typically generated. 

Figures [17] - [19] now show typical results for the good 
solvent case {q — 1) but varying the parameter qAB con- 
trolling the chemical incompatibility. The visual inspec- 
tion of the configurations (Figure [TTj) reveals little in- 
fiuence of qAB, however, and this observation is corrob- 
orated by the more quantitative analysis: the average 





FIG. 17: Snapshot pictures of bottle brush configurations for Lb = 64, q' = l,iV = 18, a = 1 and three choices of qAB, Qab = 1.0 
(a), qAB = 0.4 (b) and qAs = 0.1 (c). Monomers A, monomers B, and a straight rigid backbone are shown in black, light gray, 
and white colors, respectively. 
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FIG. 18: (a) Average number of AB pairs per monomer 
{mAB}/{Nnc) plotted vs. qab for side chain lengths N = 6, 
12, and 18. All data refer to cr = 1, = 64. (b) Specific 
heat per monoiner, Cv/{^NtIc) plotted, vs. qab for AT — 6, 12 
and 18. 
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FIG. 19: (a) Radial density distribution p(r) plotted vs. r for 
!/(, = 64, A'^ = 18, g = 1, cr = 1, / = 1 and various qab- 
(b) Normalized transverse mean square end-to-end distance 
{{Rl) + {Rl))/N^'' plotted vs. N, for Lt, = M, q = 1, a = 
1, / = 1, and three choices of qab as shown. 
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R,,„(n) 



FIG. 20: Vectors R"^{n) from the grafting site n to the xy- 
component of the center of mass of the respective chain, and 
corresponding unit vectors (denoted by arrows). For a per- 
fectly phase separated structure with the interface between A 
and B being the yz plane, for q = yl all unit vectors point 
along the negative x-axis and for a = -B all unit vectors point 
along the positive x-axis. 
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FIG. 21: Correlation function Cn measuring local "Janus 
cylinder" -type phase separation plotted versus n, for A'^ = 6 
(a), 12 (b), and 18 (c). Various choices of qab are included, as 
indicated in the figure. Curves show fits to Equation (|50|l . All 
data refer to good solvent conditions {q — 1), and Lt = 64, 
a = 1, / = 1. 
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FIG. 22: Distribution P{6) of the angle 6 between the vectors 
towards the centers of mass of subsequent (unlike) side chains 
plotted vs. e/iT, for A = 12 (a) and A = 18 (b). Various 
choices of qab are included, as indicated. 



number of A-B pairs per monomer is extremely small 
(Figure fTSk ) even for qab = 1, and hence not much 
enthalpy could be won if A-chains and B-chains avoid 
each other: due to the excluded volume interaction, very 
few nearest neighbor contacts between any non-bonded 
monomers occur in our bottle brush model. For A'^ = 18 
and QAB = 1 the total number of AB contacts per chain 
is only about 1.4, and increases with increasing N only 
very slowly. So for the range of side chain lengths ac- 
cessible in our work, no phase separation should be ex- 
pected. In the specific heat one does see a weak peak near 
0.2 < qAB < 0.4, but the height of this peak decreases 
very strongly with increasing A^. Furthermore does nei- 
ther the peak position shift with increasing A^ (as one 
would expect from Equation (|38p . if this peak would be 
a rounded precursor of the phase transition that should 
occur at Xab limit N oo) nor does the peak 

width decrease with increasing A^. Thus, it is clear that 
this peak is not an indicator of a "Janus-cylinder" -type 
phase separation in the bottle brush: rather we interpret 
it as a Schottky-type anomaly, expected from the fact 
that in our model of alternatingly grafted A-and B-chains 
at a straight line backbone with coordinates (0, 0, z) in 
the immediate environment of the backbone (e.g. at lat- 
tice sites {x,y,z) = (±1,0, z) or (0, ±l,z)) there is a 
nonzero a-priori probability of 1/4 that between the first 
monomer of the chain grafted at z and the first monomer 
of the chain grafted at 2: -I- 1 a nearest-neighbor contact 
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FIG. 23: (a) Specific iieat per monomer, Cn /{Nric), plotted 
vs. qab for A'' = 6, 12, and 18 for a = 1, Lt = 64, and 
q = 1.3087 (a Theta solvent), (b) Radial density distribution 
p{r) plotted vs. r for Lj, = 64, iV = 18, q = 1.3087, cr = 1, 
/ = 1, and three choices of qAs- 
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FIG. 24: (a) Normalized transverse mean square end-to-end 
distance (< > -I- < ii^ >)/N plotted vs. N, for L,, = 64, 
a — 1, f — 1, and q = 1.3087 (Theta solvent). Three choices 
of qAB are included, (b) Distribution P{9) of the angle be- 
tween the vectors towards the centers of mass of subsequent 
(unlike) side chains plotted vs. 0/tt, for A'^ = 18, other param- 
eters as in (a). 



occurs. The finite energy from this local contacts near 
the backbone gives rise to the peak in the specific heat. 

Also the radial density profile (Figure \Wh ) shows lit- 
tle effects of varying qab, and there is also no effect in 
the gyration radius component {Rl±) = {R^gx) + {R^gy) 
of the side chains, although a weak increase occurs in 
the corresponding component of the end-to-end distance 
(Figure [T^). In the latter figure, one can see for small 
N an even-odd oscillation, but this lattice effect clearly 
has died out for N > 10. 

In order to test for correlations measuring local phase 
separation along the backbone. Equation ([^5]) is some- 
what cumbersome to implement numerically, since for 
each z one has to find the direction of ip{z) from the 

condition that |V'(^)| is maximal (cf. Equation (|3ip ). A 
simpler and similar correlation function has been defined 
from the vectors pointing from the grafting site i to the 
center of mass of the respective side chains. Projecting 
this vector into the xy-plane and defining an unit vector 
S°' {a — A or B) along this projection (Figure (201) we 
define a correlation function C„ as follows 

C„ ^ [{S^ ■ Stn) + (5f+i • 5'f+i+J]/2 . (49) 

Here we have explicitly incorporated the alternating 
grafting ABAB... of side chains along the backbone. The 
average (. . .) in Equation (I49p includes an averaging over 



all sites {i} on which A chains are grafted, in order to im- 
prove the statistics. If perfect long range order occurs, as 
implied in Figure !^ we clearly have C„ = 1 independent 
of n, while for the case of short range order, we expect 
Cn oc exp(— n/^). Actually, considering the fact that we 
use a periodic boundary condition, we have analyzed our 
numerical data in terms of the ansatz 

C„ cx {exp(-n/0 + exp[-(Lb - n)/^]} . (50) 

Figure [m shows our data for C„ for the three choices of 
N: indeed we recognize that C„ decays to zero with in- 
creasing n, but the increase does get slower with increas- 
ing side chain length N . The scale of this correlation 
effect clearly increases with decreasing qAB- While for 
= 6 the correlation length ^ hardly depends on qABi 
for large N a slight increase of ^ with decreasing qAB is 
suggested. 

In an earlier simulatpn study of phase separation in 
binary bottle brushes, it was suggested to quantify 
the degree of separation by considering the distribution of 
the polar angle <pi from the axis to monomer Defining 
a variable af ^ which is af = 1 if monomer £ is of type 
A and zero otherwise, and similarly af , de Jong and ten 
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FIG. 25: Correlation function d { Equation (|49} } plotted 
vs. n, for Lt = 64, q = 1.3087, a ^ 1, f ^ 1, N ^ 12, 
and various choices of qab, as indicated. Case (a) refers to a 
choice of unit vectors from the z-axis to the center of mass of 
the chain that is grafted a,t z = i, while case (b) refers to a 
choice of unit vectors from the z-axis to the center of mass of 
all monomers of type a in the xy-plane ai z = i. Curves are 
fits to Equation (f50)). 



Brinke ^''i introduced a function 

P(^) = T.Y.^"^''^ + cb) (51) 

e e' 

and studied P{(t>) varying XAB- However, it turns out 
that P{4>) always has a rather complicated shape, and 
its dependence on xab is rather weak. Thus we investi- 
gated a related but somewhat different function, namely 
the histogram P{9) of the angle 9 between the vectors 
from the axis to the centers of mass of subsequent unlike 
side chains (Figure [H]). We see in this distribution P{9) 
also only a rather small dependence on xab, however: 
while for qab — 1 this distribution is essentially struc- 
tureless for tt/2 < 9 < 37r/2, for small qab = 1 a rather 
pronounced peak at « tt develops, indicating a prefer- 
ence of antiparallel orientation of subsequent side chains. 
However, we feel that such indicators as P((/') or P{9) are 
sensitive only to the presence of short range order rather 
than long range order, and hence we shall focus on the 
correlation function C„ in the following. 

Next we focus on binary bottle brushes in a Theta sol- 
vent {q = 1.3087). Figures [23] - [2S] show that the results 
are not very different from the good solvent case: the av- 
erage number of AB pairs per monomer (not shown) and 



the specific heat (Figure [23k) are hardly distinguishable 
from the good solvent case. Also the radial density pro- 
file (Figure [^5b) exhibits only minor differences, slightly 
larger densities occur inside the brush than in the good 
solvent cases, and the transverse component of the end- 
to-end distance does not exhibit much additional stretch- 
ing, when qab decreases (Figure [Mk ). The distribution 
P{9) is somewhat less flat near 9 = tt (Figure [Mb) than 
in the good solvent case (Figure [22)1 . 

Also the correlation function C„ shows again an expo- 
nential decay with n (Figure [^5k). the correlation length 
^ being rather similar to the good solvent case (Fig- 
ure I2ip . Since one might argue that our definition of 
a correlation function as given in Equation (|49p is not 
the optimal choice, and there might occur larger corre- 
lation lengths for a more clever choice of a correlation 
function, we tried a different choice which is also easy 
to compute. Namely, in the xy-plane at the index i of 
the z-coordinate we determine the center of mass of the 
monomers of type a = A or B in that plane (the choice 
of a is dictated by the type of chain grafted at z = i). 
Then we introduce an unit vector S" from the z-axis in 
the direction towards this center of mass. In terms of 
these unit vectors, which then no longer distinguish from 
which chain the monomers a in the i'th xy-plane are com- 
ing, we can again apply the definition. Equation to 
derive a correlation function. This correlation function 
is shown in Figure [^5b. One sees that the qualitative be- 
havior of both types of correlation functions is the same, 
but the decay of this second type of correlation function 
even is slightly faster than that of the correlation function 
used previously. Since normally, when one studies prob- 
lems involving a phase transition, the correlation length 
of the order parameter is the largest correlation length 
that one can find in the system, the previous definition 
(focusing on the location of the center of mass of the in- 
dividual side chains) seems preferable to us. It is possible 
of course that such differences between different ways of 
measuring correlations along the backbone remain so pro- 
nounced for small N only, where the correlation lengths 
are only of the order of a few lattice spacings, and hence 
are less universal and depend on the details of the studied 
quantity. 

We now turn to the poor solvent case (Figures [221- dS])- 
Already the snapshot picture of the bottle brush poly- 
mer conformations (Figure I26p reveals that now the side 
chains adopt much more compact configurations, due to 
the collapse transition that very-Jpng single chains would 
experience in a poor solvent. If now qab becomes 
small, one recognizes a more pronounced phase separa- 
tion along the backbone of the polymer, although there 
still is no long range order present. 

Still, neither the variation of the average number of AB 
pairs per monomer {itiab) / {Ntic) with qab (Figure [?7b) 
nor the specific heat (Figure [27b) give a hint for the oc- 
currence of a phase transition: in fact, these data still 
look like in the good solvent case. Figure [IS] Also the 
distribution P{9) and the variation of {{Rl) + (Rl))/N^^^ 



FIG. 26: Snapshot pictures of bottle brush conformations for Li, = 64, q = 1.5, A'^ = 18, o" = 1, and three choices of qAB, 
namely qAB ~ 1-5 (a), qAB = 0.6 (b) and qAB = 0.1 (c). Monomers A, monomers B, and a straight rigid backbone are shown 
in black, light gray, and white colors, respectively. 



with N and qab are qualitatively similar to the results 
found for good solvents and Theta solvents, and there- 
fore are not shown here. More interesting is the monomer 
density profile (Figure [28]) . While for = 6 it still has 
the same character as in the previous cases, for = 12 
and = 18 we recognize an inflection point. In fact, 
for a collapsed polymer bottle brush we expect a profile 
exhibiting a flat interior region at melt densities {p{r) 
near to p = 1), then an interfacial region where p{r) 
rapidly decreases towards zero; in the ideal case reached 
for q — > oo, the profile should even be a Heaviside step 
function, in the continuum limit, p{r) — 9{h — r). Obvi- 
ously, for g = 1.5 we are still far from this behavior even 
for N — 18, but we can identify an interface location at 
about Tint ~ 3.5, and in the interior of the cylindrical bot- 
tle brush (for r < Tint) the monomer density corresponds 
at least to a concentrated polymer solution. 

Figure then shows again plots of C„ vs. n, and the 
resulting fits to Equation (|50p . The behavior is qualita- 
tively similar to the previous cases of good and Theta 
solvents again, but now the decay of C„ with n is clearly 
much slower, indicating a distinctly larger correlation 
length. This result corroborates the qualitative observa- 
tion made already on the basis of the snapshot pictures. 
Figure \26\ and quantifies it. 

Our findings on the local character of phase separa- 
tion in binary bottle brush polymers are now summa- 
rized in Figure [30| where plots of the inverse correlation 
length 1/^ versus the inverse Flory-Huggins parameter 
Zc/x = kBT/{eAB - e-AA) = l/\n[q/qAB\ are shown (cf. 
Equation ([28]) and remember our choice caa — ^bb)- 
This choice of variables is motivated by the analogy with 
the XY model in cylindrical geometry. Equation (|45p ; on 
the basis of this analogy, we would expect that cx T, 



i.e. we expect straight lines that extrapolate through 
zero as T ^ 0. |— ■ 

Figure [3D] demonstrates that this analogy clearly 
is not perfect. Rather the data are compatible with an 
extrapolation 

r^{N,T)^C^{N,q,Q)+CN{q)T , T->0, (52) 

where C^iq) is a coefficient that decreases with increas- 
ing TV and depends on solvent quality and with a nonzero 
intercept ^~^(N,q,0) implying that even in the ground 
state (T — 0) there is lack of long range order! 

A very interesting question concerns the dependence 
of ^ {N, T) on chain length. Figure [ST] presents a plot 
of ^~^(iV, (J, 0), as estimated from extrapolation of the 
data shown in Figure [30] to qAB = 0, as a function 
of 1/A^. Indeed the data are compatible with a rela- 
tion ^~^(A^, g,0) oc \/N, implying long range order for 
iV — > oo in the ground state. In principle, carrying out 
the same extrapolation at nonzero temperature and test- 
ing for which temperature range £,-'^{N oo,T) starts 
to be nonzero, one could obtain an estimate for the criti- 
cal point of intramolecular phase separation, that should 
occur (and then is well defined) in the limit of infinite 
side chain length, N ^ oo. Unfortunately the accuracy 
of our estimates for ^{N,T) does not warrant such an 
analysis. 

The reason for this surprising persistence of disorder 
at low temperatures is configurational entropy, of course. 
As long as g < oo, the bottle brush polymer is not fully 
compact, and when the density is small enough inside 
the bottle brush (as it still is in the cases shown in Fig- 
ure [25]), the two types of chains A, B can avoid mak- 
ing binary contacts (as is evident from Figures [TH] and 
[?7l) : {rriAB) — > as qAB 0. Hence there is a per- 
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FIG. 27: (a) Average number of AB pairs per monomers, 
< ruAB > /(Nric), plotted vs. qab for side chain lengths 
iV = 6, 12, and 18. All data refer to a = 1.5, Lt = 64. (b) 
Specific heat per monomer, Cv / [Nric), plotted vs. qas for 
iV = 6, 12, and 18. 




FIG. 28: Radial density profile p(r) plotted vs. r for Lb — 64, 
qAB = 0.1, q — 1.5, cr = 1, / = 1, and three choices of A'' = 6, 
12, and 18. 



feet avoidanee of energetieally unfavorable contacts in 
the ground state, and the finite correlation length ^{N, 0) 
and nonzero entropy of the groundrState are not a conse- 
quence of "frustrated interactions^UnpLas in spin glasses, 
random field spin models, etc. 11221. liOal the present 
problem, if both q and N are not too large, a ground 
state (for q ^ 0) occurs where the structure of the bottle 
brush is NOT a cylinder with an interface separating an 
A-rich and a B-rich domain, as hypothesized in Figure[T3] 
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FIG. 29: Same as Figure [21] but for poor solvent conditions 
{q = 1.5); for n = 6 (a), n = 12 (b), and n = 18 (c). 



Thus a (coarse-grained) cross section of the bottle brush 
in the xy-plane is NOT a circle, separating A-rich and 
B-rich regions by a straight line, but rather looks like the 
number 8, i.e. a dumbbell-like shape, where a (more or 
less circular and more or less compact A-rich region) oc- 
curs for X < 0, a similar B-rich region occurs for a; > 0, 
but no monomers (apart from the backbone monomer) 
occur at the y-axis at a; = 0. The orientation of the 
X-axis can fluctuate as one moves along the z-axis, and 
unlike Figure [15] (twisting an interface in a cylinder along 
the z-axis clearly involves an energy cost, and this is de- 
scribed by the helicity modulus in the XY-model analogy, 
Equation ([iS]) ) the free energy cost of this structural dis- 
tortion is outweighed by the configurational entropy gain, 
at least for q and N not too large. It is an unresolved 
question whether some critical value qdN) exists, where 
^~^{N,q,0) vanishes, and a long range ordered ground 
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FIG. 31: Inverse correlation length plotted vs. 1/A^, in the 
limit qAB — > 0, for poor solvent conditions {q = 1.5), and 
three choices of the backbone length Lj,, to exclude that the 
data are strongly affected by finite size artifacts. The straight 
line indicates that the data are compatible with a simple linear 
relation, ^"^(iV,0) oc 1/iV. 
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FIG. 30: Plot of the inverse correlation length 1/^ versus 
the inverse Flory Huggins parameters, Zc/Xi for good solvent 
conditions {q — 1), case (a), Theta solvent conditions {q — 
1.3087), case (b), and poor solvent conditions (g — 1.5), case 
(c). Data for three chain lengths A'^ = 6, 12, and 18 are 
included throughout. In several cases two or three choices of 
backbone lengths Lt — 32, 48 and 64 are included as well. All 
data refer to the choice of one grafted polymer per grafting 
site (cr = 1). 



state occurs. Another unresolved question is, whether 
(another?) critical value q'dN) exists, where the charac- 
ter of the ground state changes such that the local cross 
section of the binary bottle brush changes from an 8- 
shaped to a circular density distribution. Actually, the 
problem of alternatingly grafting A-chains and B-chains 
along a line is equivalent to grafting symmetric AB di- 
block copolymers along a line, such that the junction 
points of the diblocks form a straight line. This prob- 
lem is the lower-dimensional analog of diblocks grafted 
with their junction points to the flat interface between 



an unmixed binary (A,B) homopolymer bleipd^ for this 
problem it is well-known (see Werner et al. '^231 for ref- 
erences) that the shape of the diblocks is dumbbell like, 
the A-block being stretched away from the interface, in 
order to be embedded in the A-rich phase underneath of 
the interface, and the B-rich block also being stretched 
away from the interface, in order to be embedded in the 
B-rich phase above the interface. 

As a consequence of our discussion, we call into 
question tli©, idea^of the "Janus cylinder" -type phase 
separation ISa. ISIl ISSl and propose as an alternative pos- 
sibility (Figure [5^ the "double cylinder" (with cross sec- 
tion resembling the number 8). Which of these cross- 
sectional structures occur will depend on the interac- 
tion parameters eAA — ^bb, and cab, of course. If 
the strength of the attractive interactions \eAA\ = \f-BB\ 
exceeds the strength of the repulsion \eAB\ sufficiently 
much, it is clear that for T ^ the Janus cylinder type 
phase separation will occur, while in the opposite limit, 
when I CAB I exceeds \€aa\ — \^bb\ sufficiently much, the 
double cylinder geometry will win. In the macroscopic 
continuum limit a comparison of the respective surface 
energies would imply that the double cylinder geometry, 
which avoids an AB interface, but requires for the same 
volume taken by A and B monomers a surface area of 
pure A and B that is larger by a factor of ^/2 than for 
the "Janus cylinder" (Figure [5^ . becomes energetically 
preferable for 

iV2 - l)TT\eAA\ = l.30l\eAA\ (53) 

Since in the present work we consider the limit qab = 
exp(— EAs/fesT) ^ at fixed q = exp(— eAA/fcsT"), it 
is clear that in our case Equation ([55]) is fulfilled. How- 
ever, for the "most symmetric" choice of interaction 
parameters, e^s — —^aa — ^^bb, the conclusion would 
be different. For lattice models Equation ([55)1 is ques- 
tionable since very compact configurations of collapsed 
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polymers must respect the lattice structure, and different 
geometrical factors, depending on the type of the lattice, 
in the inequality Equation ([53|l may occur. 



VII. Conclusions and Outlook 

In this article, we have restricted attention ex- 
clusively to static conformational properties of very 
long bottle brush polymers with a rigid back- 
bone. We feel that this "simple" limiting case 
needs to be understood first, before the very inter- 
esting extensioa, to. Ae— ca«€|-Qjf|^aidble. psemifiexible 
backbones, ^ . y M |3i|, y , Isi Isi Isi Isi y, y, Isi y or the 
question of the crossoveipbetween a bottle brush poly- 
mer and a star polymer, can be correctly addressed. 
Thus the latter two problems were not at all considered 
in this paper, and hence w9-aJs©,(i©,not wish to comment 
on the recent controversy li^ lla. Ua concerning the cor- 
rect interpretation of experiments on the overall linear 
dimensions of bottle brush polymers. Thus, the focus of 
the present article, as far as one-component bottle brush 
polymers are concerned, is the conformation of the side 
chains. We recall that this information is experimentally 
accessible, if one prepares bottle brushes with a single 
arm being deuterated, while all remaining arms of the 
polymer remain protonated, to allow a study of the static 
structure factor S{q) of a single arm by elastic coherent 
neutron scattering. It is clear that such experiments are 
very difficult, and hence no such experiment of this type 
is known to the present authors yet, but clearly it would 
be highly desirable to obtain such experimental informa- 
tion. 

As discussed in the first part of the paper, one can find 
in the literature rather diverse concepts about the con- 
formations of the side chains of a bottle brush polymer 
under good solvent conditions and high grafting density. 
One concept assumes that the cylindrical volume that 
the bottle brush occupies can be partitioned into disks, 
such that each disk contains just one polymer chain con- 
fined into it, no other chain participating in the same 
disk. A simple geometric consideration, reviewed in the 
first part of the present paper, then yields the following 
predictions for the linear dimensions of the chain, as a 
function of grafting density cr and side chain length N 

Rg^ oc a^/^N""/*, Rgy cx a^/^N^/\ Eg, oc a'^ . (54) 

Remember that we choose the z-axis as the direction of 
the rigid backbone, the x-axis is oriented normal to the 
z-axis towards the center of mass of the chain, and the y- 
axis is perpendicular to both x- and z-axes. Thus the 
chain conformation has a quasi-two dimensional char- 
acter, and there is no stronger stretching of a chain 
in radial direction than in tangential direction, since 
Rgy/Rgx = coust (indcpendcnt of both N and cr). 

Clearly, this quasi-two-dimensional picture is not very 
plausible, and more popular is an extension of the Daoud- 



Cotton blob picture for star polymers to the present 
case. While for a star polymer the blob radius simply 
increases proportional to the distance from the center, r, 
and there is no geometrical difficulty to densely pack the 
conical compartments resulting from dividing a sphere 
into / equal sectors, such that each sector contains a sin- 
gle arm of the star, with a sequence of spheres of increas- 
ing size, it has been argued in the literature that for cylin- 
drical geometry the blob radius scales as ^(r) cx (r/a)^/^. 
While in the literature no mentioning of a non-spherical 
blob shape (or, equivalently anisotropic local screening of 
the excluded volume interaction) is found, we have em- 
phasized here the geometrically obvious fact that a dense 
filling of space with blobs in the cylindrical geometry ap- 
propriate for a bottle brush requires that the blobs have 
the shape of ellipsoids with three different axes, propor- 
tional to r, (r/cr)^/^, and cr^^, respectively. While this 
picture does not alter the prediction for the stretching of 
the chain in the radial x-direction found in the literature, 
Rg,^ cx cr(i-'')/(i+'')iV2'^/(i+'') « cri/47v3/4^ different pre- 
dictions result for the other two linear dimensions. While 
the regular Daoud-Cotton picture would yield for Rgy 

the size of the last blob, Rgy — f (r ~ h) = [h/aY^'^ cx 
^-y/(i+v)j^u/{i+u) ^ ^-3/8^3/8^ gjj^ instead 

Rgy CX ^-(2>'-l)/(l+-)^3./(2.+2) ^ ^-1/16^9/16 ^ (55) 
Rg,_ cx ^-(2--l)/(l+-)iV-/(2-+2) « ^5/8^3/8 _ (gg) 

Equations (|55p and ([55]) make use of both the anisotropic 
blob linear dimensions and random walk-type arguments 
(cf. Equations (|13p - ([28l) ) and are new results, to the best 
of our knowledge. 

Also the crossover scaling towards mushroom behav- 
ior has been considered. Equations (P0)l - (|23p . and Monte 
Carlo evidence for this crossover scaling description was 
obtained, from extensive work using the PERM algo- 
rithm. Although rather large bottle brushes were simu- 
lated (backbone length up to Li, — 128, avoiding free 
end effects by periodic boundary conditions, side chain 
length up to = 2000 for grafting density cr = 1 , see Fig- 
ure[2]), it was not possible to reach the asymptotic regime 
of strong side chain stretching where Equations (j55p and 
(|56p hold, however. One reason why for moderate graft- 
ing densities even fairly long side chains can avoid each 
other with small significant amount of stretching is the 
fact that the natural shape of a self avoiding walk con- 
figuration is an elongated ellipsoid, with three rather dif- 
ferent eigenvalues of the gyration tensoij^-So side chains 
can to a large extent avoid each otherl^S by orienting 
themselves such that the longest axis of the ellipsoid is 
oriented in the radial x-direction and the smallest axis 
in the backbone z-direction. This implies that very large 
values of aaN'^ are needed to obtain significant stretch- 
ing. This difficulty to verify any such scaling laws in our 
simulations, which were only able to explore the onset of 
stretching away from simple mushroom-type behavior of 
the side chains and not the strongly stretched behavior. 
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FIG. 32: Possible ground structures of binary bottle brush polymers, showing the cross section perpendicular to the backbone 
(dot) in the continuum limit, for the symmetric case eAA = ess- Left part shows a double cylinder, each cylinder has a radius 
h. Right part shows a Janus cylinder (with the same volume this requires a radius h' — y/2h. The surface energy cost of the 
double cylinder is AhirLb \ €aa |, while the surface energy cost of the Janus cylinder is 2y/2hLb{'!v \ eAA \ + | tAB [)■ (b) Possible 
generalizations of the double cylinder structure to the asymmetric case eAA ess- The left part assumes ess = 0, so the B 
chains are in good solvent, only the A-chains collapse. The right part assumes | ess |>i eAA |, so that the collapsed B-rich 
cylinder is denser than the A-rich cylinder. 



suggests that no such scaling behavior should be observ- 
able in the experiments as well: given the empirical fact 
that one bond of a coarse-grained lattice mpclel--porre- 
sponds to n = 3 — 5 chemical monomers, l^'lSal it is 
clear that the experimentally accessible side chain lengths 
do not exceed those available in our simulation. Thus, 
scaling theories of bottle brushes are unfortunately of 
very restricted usefulness for the interpretation of either 
experiments or simulations. The only firm conclusion 
about the theories mentioned above that we like to make 
is that there is no evidence whatsoever for the quasi- 
two-dimensional picture, since we do see a decrease of 
Rgy/Rgx wlth iucrcasiug N. 

Turning to the problem of intramolecular phase sepa- 
ration in binary (AB) bottle brush polymers, we have ex- 
amined the proposal jthcjt-a "Janus cylinder" -type phase 
separation occurs. '^'^■'^ This idea is questionable for 
several reasons (i) In a quasi-one-dimensional system, no 
sharp phase transition to a state with true long range 
order can occur; at most one can see a smooth increase 
in the corresponding correlation length, as the temper- 
ature is lowered (or the incompatibility between A and 
B is enhanced, respectively). Thus we have defined suit- 
able correlation functions and studied the variation of the 
corresponding correlation length as function of the chain 
length N and the parameter q^B = exp[— e^s/fc^T] con- 
trolling the incompatibility {^ab is the repulsive energy 
encountered when two neighboring lattice sites are oc- 
cupied by monomers of different kind). Three choices 
of solvent quality (taken symmetric for both A and 
B through the choice caa = ^bb) were considered, 
q — eyip{~eAA/kBT) — 1 (good solvent), q — 1.3087 
(Theta solvent) and q = 1.5 (poor solvent). It was found 
that in all cases the correlation length suited to detect 
Janus-cylinder type ordering increases rather weakly as 
qAB — > 0, approaching finite values even for qAB = 0, for 



finite side chain length N, while for TV ^ oo an infinite 
correlation length is compatible with the data. No evi- 
dence for the predicted critical points and their scal- 
ing behavior with N (Equations [32ll38| could be detected, 
however. The gradual establishment of a local phase sep- 
aration lacking a sharp transition is conmatible with ob- 
servations from a previous simulation, however, (ii) 
Depending on the relation between the energy parame- 
ters eAA and eAB (cf. Equation the "local" phase 
separation (considering a slice of suitable thickness per- 
pendicular to the backbone of the bottle brush to obtain 
suitable coarse-grained densities pa,Pb on mesoscopic 
scales) for temperatures T — > may have two differ- 
ent characters (Figure [5^): the "Janus cylinder", which 
contains an interface between A-rich and B-rich phases, 
competes with the "double cylinder" . In the latter struc- 
ture, there is no extended AB interface, A-chains and B- 
chains meet only in the immediate vicinity of the back- 
bone where they are grafted. Thus, the coarse-grained 
density distribution under poor solvent conditions in a 
slice has the shape of the number 8, where in the upper 
part of the 8 we have the B-rich phase and in the lower 
part we have the A-rich phase. Of course, when we ori- 
ent the X-axis (Figure such that it is parallel to the 
vector connecting the center of mass of the A-rich region 
to the center of mass of the B-rich region in the slice, the 
orientation of the x-axis for both structures in Figure [32] 
can randomly rotate when we move along the backbone 
(z-direction) due to long wavelength fluctuations, cf. Fig- 
ure [151 and hence the comments about finite correlation 
lengths of this phase separation apply here as well. Ac- 
tually, the numerical results of the present Monte Carlo 
simulations give clear evidence that local phase separa- 
tion of "double cylinder" -type rather than "Janus cylin- 
der" type is observed, since the number of AB-contacts 
tends to zero. 
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FIG. 33: Correlation function Cz (Equation JST])) plotted 
vs. z, for Lb = 64, q = 1.5, a = 1, / = 1, iV = 6 (a) and 
A'' = 18 (b), and various choices of qab, as indicated. 



Thus, varying the solvent quality and chemical incom- 
patibility, one can influence the character of local in- 
tramolecular phase separation, and one can control the 
correlation length over which the vector characterizing ei- 
ther interface orientation or the axis of the local dumbbell 
is oriented in the same way along the z-axis. Of course, 
in real systems one must expect that the solvent quality 
for the two types of chains will differ (Figure [32b). Then 
asymmetric 8-shaped local structures will result: e.g., if 
the solvent is a good solvent for B but a poor one for A, 
but the incompatibility between A and B is very high, we 
expect a structure where the A-chains are collapsed in a 
cylinder with the backbone on the cylinder surface, and 
from there the B-chains extend into the solution like a 
"flower" in the cross section (Figure [5^. left part). Con- 
versely, it may happen that the solvent quality is poor 
for both A and B, but nevertheless different, so the den- 
sities of both cylinders and hence their radii will differ 
(Figure [5^. right part). Similar asymmetries are also of 
interest if the local character of the phase separation is 
of the Janus cylinder type: then in the cross section the 
A-rich and B-rich regions will take unequal areas, rather 
than equal areas as shown in Figure [32k . and the inter- 
face will be bent rather than straight. Similar effects will 
occur when the chain lengths Na,Nb differ from each 
other, or the flexibilities of both types of chains are dif- 
ferent, etc. 



An aspect which has not been discussed here but which 
is important for the binary bottle brush under poor sol- 
vent conditions is the question whether or not the density 
in the collapsed bottle brush is homogeneous along the 
z-direction (so one really can speak about "Janus cylin- 
ders" or "double cylinders" , Figure [5^ , or whether it 
is inhomogeneous so the state of the system rather is a 
chain of "pinned clusters" . For the one-component bptr 
tle brush, this question was discussed by Sheiko et al.,'Sfill 
from a scaling point of view, but (unlike the related prob- 
lem of bjraah-eiuster|-tKinsition of planar brushes in good 
solvents^ii^ [^^^ [^^^ [^^^ ) we are not aware of any simula- 
tion studies of this problem yet. In order to understand 
the structures of the binary bottle brush in the poor sol- 
vent, we count the number of monomers M(z) (irrespec- 
tive of A or B) in each xy-plane, and calculate the nor- 
malized correlation function Cz of dM{z) — M(z) — N 
along the backbone, 

Cz =< SM{i + z)5M{i) > / < 5M{if > . (57) 

Here the average < . . . > includes an averaging over all 
sites {i} on which chains are grafted, in order to improve 
the statistics. In Figure [33l we see that this correlation 
actually is negative for 2 < ^ < 8, indicating some ten- 
dency to pinned cluster formation. 

Before one can try to understand the real two- 
component [—bottle brush polymers studied in the 
laboratory, 111^11131 two more complications need to be 
considered as well: (i) the flexibility of the backbone; (ii) 
random rather than regular grafting along the sequence. 
While the effects due to (i) were already addressed to 
some extent in an earlier simulation, randomness of 
the grafting sequence has not been explored at all. HjOWi- 
ever, in the study of binary brushes on flat substrates^ii^ 
it has been found that randomness in the grafting sites 
destroys the long range order of the micij©-«hase sepa- 
rated structure, that is predicted to occur^ii^ for a per- 
fectly periodic arrangement of grafting sites. In the one- 
dimensional case, we expect the effects of quenched disor- 
der in the grafting sites to be even more dominant than in 
these two-dimensional mixed polymer brushes. Hence, it 
is clear that the explanation of the structure of one- and 
two-component bottle brush polymers still is far from be- 
ing complete, and it is hoped that the present article will 
motivate further research on this topic, from the point of 
view of theory, simulation and experiment. 
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